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This thesis presents an interpretation of the Theory of 
Biorhythms and develops a combined statistical model for a 
statistical validation of the theory. The model is then used 
to investigate a set of intellectual data and a set of acci- 
dent data. The 'intellectual data is based upon 112 academic 
grades taken from four postgraduate school students over a 
14 month period. The accident data, taken from insurance 
abstracts, consists of 66 accidents which occurred at a pulp 
plant in British Columbia, Canada. 
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I. INTRODUCTION 



It is currently suggested in scientific theory that life 
functions of all living organisms are characterized by some 
form of cyclic nature. Many of the cyclic rhythm theories 
based on observed human performaiice have been accepted as 
fact while others, lacking strong scientific and mathematical 
proof, have been treated with some skepticism. 

The most well-known and best documented of observed 
human isochronic phenomena are the circadian or daily rhythms. 
Circadian rhythms vary in length from twenty to twenty- eight 
hours depending upon the individual and environmental cir- 
cumstances. Included in these daily or circadian rhythms 
are life function or regulatory cycles, biochemical cycles, 
and socio-cultural , work-eat-sleep cycles. Extensive re- 
search in the area of circadian rhythms, as reported by 
Squibb Pharmaceutical Co. [1972], concluded that regulatory 
and biochemical cycles are extremely rigid. Lengthy experi- 
mentation with circadian work-eat-sleep cycles has been con- 
ducted by a variety of individuals and institutions including 
the Max Planck Institute in Germany [Palmer] . Their research 
originated from two contesting theories. One theory states 
that the circadian patterned work-cat-sleep cycle is a prod- 
uct of socio-cultural customs while the other theory suggests 
it resulted from the interaction of the multitude of circadian 
biochemical and regulatory cycles [Luce, 1973]. Scientists 
have been conducting experiments in this area for the past 
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thirty- four years. The experimental prerequisite to remove 
artifical and natural indicators of time, such as clocks 
and the rising and setting of the sun, stimulated Dr. 

Nathaniel Kleitman to conduct his original cave studies, in 
Mammouth Cave, Kentucky [Luce, 1972c] . In 1962, further 
studies were conducted by Michel Siffre who spent sixty- 
three days on a subterranean glacier in Scarasson Cavern in 
the French Alps. Finally, in 1966, Dr. John Mills analyzed 
David Lafferty's one-hundred and twenty-seven day isolation 
in a cave in Cheddar, England. The most controlled experi- 
ments were conducted at the Max Planck Institute where 
eighty-five subjects lived for three weeks in luxurious 
underground efficiency apartments. 

All research yielded the same results that although man 
can artifically control the length of his work-eat-sleep 
cycle, the regulatory and biochemical cycles maintain their 
ritualistic twenty-four hour schedules. Research has also 
indicated that the regulatory and biochemical circadian rhythms 
influence physiological responses to environmental stimuli. 

For example, the body reacts differently to such foreign 
agents as noise, antibiotics, hypnotics, anaesthetics, and 
liquor at different times of the day [Squibb, 1972]. It has 
often been suggested that circadian rhythms are a partial 
explanation for the commonly observed circadian post-surgical 
death cycles and spontaneous birth cycles. Since circadian 
rhythms are constantly observed, the fact that they do exist 
is generally accepted, but, in some cases, why they exist can 
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only be hypothesized since the precise roles of many body 
functions remain a mystery to the medical community. 

Cycle analysis has not been limited to the study of 
. circadian rhythms. Recently behavioral scientists have be- 
come involved in the analysis and interpretation of patterns 
of overt human performance. A considerable amount of work 
has been conducted by physiologists and psychologists in 
the area of cyclic phenomena of human behavior. Development 
and application of several emotional measurement devices has 
led to various mood cycle discoveries. Tedious experimen- 
tation and research by Dr. Max Liischer of Germany led to the 
original development of the Liischer Color Test in 1947 which 
has since been extensively refined. Its validity is based 
upon several different statistical analyses [Scott, 1971]. 
Adaptability, ease of application, and lack of rigidity has 
made it an ideal and instrumental tool of many behavioral 
scientists, from psychologists to criminologists. Although 
investigation and clinical observations by Drs. Halberg and 
Panofsky revealed a wide spectrum of mood cycles, they re- 
ported the most statistically significant periodicity was 
seven days [Luce, 1972b, Brinker, 1973]. Further mood cycle 
study was conducted by J. A. Dorland [Brinker, 1973]. Using 
the Dorland Mood Rating Scale, a categorical, subjective, 
emotional, measurement device, Dorland also noted the pres- 
ence of a strong seven day mood cycle. Additional investiga- 
tion of emotional cycles has led to various other mood scales 
such as Hersey's Rating Scale for Emotions [Brinker, 1972]. 
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Additional research in the areas of emotional distur- 
bances led to more bizarre theories. Statistical studies, 
conducted by Dr. Arnold Lieber of the University of Miami's 
Psychiatry Department, based on 1800 murders in Dade County, 
Florida, show a cyclic murder rate coincidental with the 
rise of the full moon [Hackler, 1969]. Doctors from the 
Douglas Hospital in Montreal, Canada have discovered cor- 
relation between cycles of intense hostility and violence 
among mental patients and geomagnetic disturbances caused 
by sunspots. Many physicians such as Dr. M. Arborlieus of 
Halmstead, Sweden, believed that the awareness, proper in- 
terpretation, and application of physiological and psycho- 
logical behavior patterns can be of significant importance 
in accurately diagnosing certain diseases such as ulcers, 
cancer, and psychosomatic illness [Luce, 1972a]. 

One of the older cyclical theories, based on observation 
of human behavior, is the theory of biorhythms developed in 
the late 19th century by Dr. Herman Swoboda of Switzerland 
and Dr. Wilhelm Fliess of Germany [Thommen, 1973]. Based 
on research in psychology and periodicity, the theory states 
that man is characterized by three inherent cycles called 
biological rhythms or biorhythms. The three cycles, which 
include a 23-day physical biorhythm, a 28-day emotional 
biorhythm, and a 33-day intellectual biorhythm, are indica- 
tive of the potential performance capacity of an individual. 
It is hypothesized that all three biorhythms begin at birth 
and terminate with death. Swoboda and Fliess' cyclic 
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hypotheses were based on pure observation. Since the 
original work of Swoboda and Fliess, various people have 
explored biorhythms resulting in a variety of conclusions 
and interpretations. 

Wallenstein and Roberts [1973] tested biorhythm validity 
against the outcome of the 1972 UCLA-USC football game along 
with selected historical events including Custer's prognosis 
at Little Big Horn and Babe Ruth's sixty homerun season. 
Thommen [1973] attacked validity in a similar fashion by 
analyzing specific historical events biorhythmically fol- 
lowed by. a subjective interpretation of the results. Willis 
[1972] analyzed data which included hospital records, single 
car accident fatalities and athletic events. He also dis- 
cussed several cases of current application including the 
supposed success of the Ohmi Railway Company of Japan in 
achieving 2,000,000 accident free kilometers since the im- 
plementation of biorhythms in driver management. Although 
the chosen specific examples have corresponded favorably to 
biorhythmic prediction, this type of analysis does not con- 
stitute a proper mathematical validation of the theory. 
Closer scrutiny of Ohmi's accomplishments may not be so sur- 
prising if an analyst properly applied stochastic processes. 
A time or spatial poisson counting process , which counts 
events through time in a poisson fashion, could have been 
applied to predict the probability of achieving 2,000,000 
accident free kilometers. Using an estimated accident rate 
based on past performances it may be entirely possible that 
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this could have occurred with relatively high probability, 
prior to the use of biorhythms. The following brief mathe- 
matical formulation illustrates the application of a sto- 
chastic process to the investigation of the Ohmi situation. 
To be complete it requires the appropriate data, from the 
Nagahama Service Center, needed to estimate the accident 
rate. However, 

P(n accidents occur in the interval (s,s+t) = 



P (N (s+t) - N(s) = n) = P (N (t) = n) = 



e" Xt (Xt) n 

n! 



For the Ohmi situation: n = 0 



t = 1.083 years 
X = ft accidents/given 



P (N(l . 083) = 0) = 



e (Xt) 



-1.083X 

e 



time period 



Details of calculations, theoretical assumptions, and 
the associated postulates are contained in "Introduction to 
Probability Models" by Ross [1973]. 

The absence of a methodological approach, justification 
for chosing data categories and acceptable mathematical tests 
resulting in statistically significant quantifiable results, 
has led to incomplete results and subjectively interpreted 
conclusions . 
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II. STATEMENT OF PROBLEM 



The primary objective of the thesis was to establish 
the basis for a statistical validation of the Theory of 
Biorhythms. It can be assumed that the general lack of 
acceptance and understanding of the theory is the cause of 
the -difficulty to locate literature on the subject. The 
literature available is scarce and the subject matter is 
treated subjectively; consequently, a statistical analysis, 
in depth, would require careful interpretation of the theory. 
The theory has not been explicitly stated in detail in any 
of the literature while conflicting assumptions supporting 
each version of the theory depend upon the author. One 
problem area has been determining what a biorhythm actually 
measures. Willis' [1972] interpretation suggests that bio- 
rhythms are indicators of an individual's relative degree of 
potential performance capacity for a given time frame. Other 
potential problem areas include the interpretation of dif- 
ferent biorhythm interactions when analyzing observed phe- 
nomena and proper classification of events as either physical, 
emotional, or intellectual. For example, Thommen [1973] 
classifies agility as an intellectual event while Willis im- 
plies that it is a physical event. Improper classification 
results in erroneous conclusions. It is often the case that 
emotional interactions contribute heavily to outcomes of 
athletic contests, consequently, athletic contests can not 
be analyzed strictly as a physical event. Most of the 
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literature available consists of short interpretations of 
biorhythms and attempts to draw conclusions about the theory 
based upon subjective correlations between biorhythmic pre- 
dictions and special cases of actual observed phenomena. 

The lack of a strong statistical foundation is one major 
stumbling block to universal acceptance of the Theory of 
Biorhythms. Statistical validation of the theory could pos- 
sibly lead to revolutionary achievements in education, 
athletic strategies, and industrial safety. 
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III. INTERPRETATION OF THE THEORY OF BIORHYTHMS 



Human behavior is characterized by three inherent cycles 
called biological rhythms (biorhythms) which include a 23- 
day physical biorhythm, 28-day emotional biorhythm, and a 
33-day intellectual biorhythm. Each biorhythm represents 
the individual's relative degree of potential performance 
capacity. It neither predicts what will happen in the 
future nor what happens in the present. Rather, biorhythms 
are indicators of the individual's physical, emotional, and 
intellectual performance capacity. In theory, it measures 
whether the potential capacity is greater than or less than 
some neutral standard on particular days within the cycle. 
That is, biorhythms measure the day to day variation in 
performance potential. The neutral standard represents the 
individual's neutral or average performance potential and 
is different for different individuals. A realization of 
performance potential occurs when a subject encounters a 
challenging situation. A challenge allows him to identify 
his full potential. Average days are days on which the in- 
dividual does not encounter challenging situations which 
require a display of emotional, intellectual, and/or physical 
potential. For this reason, the individual should be aware 
of the different catalysts that may trigger potential re- 
actions capable of being measures by a biorhythm. 

Figure 1 is a graphical interpretation of The Theory of 



Biorhythms . 



Physical Biorhythm 




Emotional Biorhythm 




Intellectual Biorhythm 



C 




33 



C = Critical Days = A day on which the individual's potential 
capacity is considered unstable. On critical days, it is 
hypothesized that individuals can either have an extraor- 
dinary increase in performance or an extraordinary decrease 
in performance. 

1 = Good Performance Days = Days on which there is a greater 
potential for above average performance. These days occur 
within the first h cycle. 

II = Poor Performance Days - Days on which there is a greater 

potential for below average performance. These days occur 
within the second % cycle. 

Average Performance Day = A day on which the individual has 
not encountered a challenging situation. Average days 
can occur anywhere within the cycle. 

Figure 1. Graphical Illustration of .Biorhythms. 



( 
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Thornmen [197 3] suggests the following subcategories are 
contained in each cycle. 



Physical Cycle 

physical strength 

endurance 

energy 

resistance 

confidence 



Emotional Cycle 

sensibility 

nerves 

feelings 

intuition 

cheerfulness 

moodiness 

creative ability 



Intellectual Cycle 

intelligence 

memory 

mental alertness 
logic 

reasoning power 
reaction 
agility 
ambition 



Table I. Subcategories of Each Biorhythmic Cycle. 



IV. SELECTION OF THE STATISTICAL MODEL 



Although a variety of the statistical techniques avail- 
able can be used to investigate biorhythm theory, model se- 
lection is dictated by characteristics of the data, as well 
as the appropriateness and quality of the model. 

A. CHARACTERISTICS OF THE DATA 

Three approaches are available. The first approach 
consists of collecting data resulting from a phenomena which 
is predictable from biorhythm theory and then selecting a 
statistical technique whose assumptions are satisfied by 
the data and whose results most appropriately illustrate the 
validity of the theory of biorhythms. In using this approach 
it should be obvious that the choice of the technique is in- 
fluenced by the characteristics of the data, which must con- 
form to the assumptions of the statistical model. For 
example, the most common assumption of parametric tests is 
homogeneity of sampling units. Although it may be possible 
to restore homogeneity through transgeneration, the use of 
non-homogenous data in a parametric test leads to invalid 
results. Consequently, the absence of proper data, which 
can force the use of a less appropriate statistical test, 
may result in inconclusive evidence supporting the theory. 

The second approach consists of selecting the most ap- 
propriate statistical tests which best illustrate the validity 
of the biorhythm theory, then collect data conforming to the 
assumptions of that test. 
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A third approach was considered out of necessity. Data 
was often difficult to obtain which both illustrated par- 
ticular aspects of the theory and conformed to the assump- 
tions of the chosen statistical model. Consequently, it 
was necessary to transgenerate raw data into categorical 
variables and then proceed with a non-parametric test. 

B. APPROPRIATENESS OF STATISTICAL MODEL 

Appropriateness of the statistical test is determined 
by the objective of the thesis. Various parametric and 
non-parametric tests were considered. The robust t-test 
analysis of the variance, and their non-parametric counter- 
parts would have been inappropriate since they allowed anal- 
ysis of only specific segments of the cycle. Consequently, 
even if the results were favorable, the evidence would be 
inconclusive. For example, depending upon homogeneity 
qualities of the data, analysis of the variance could have 
been used to compare the statistical differences between 
biorhy thmical ly predictable critical days and the actually 
observed critical days. However, this technique would not 
allow a simultaneous investigation of the statistical sig- 
nificance of the occurrence of other categorical events be- 
tween the biorhythmic critical days. The results of the 
statistical models must completely support the objective 
of the thesis. Therefore, tests which offer only a partial 
analysis were considered inappropriate. 

An appropriate statistical model for testing the validity 
of biorhythms should test the existence of a relationship or 
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dependency between observed good performance days and bio- 
rhythmically predictable good performance days and observed 
poor performance days and biorhythmically predictable poor 
performance days. Figure 2 illustrates this relationship. 

It is intuitively appealling that a relationship exists 
since most of the good performance days occur during the 
first % cycle and most poor performance days. occur during 
the second h cycle. An appropriate statistical model should 
be able to determine the presence of a relationship or de- 
pendency between theoretically predictable days and observed 
days as illustrated. If the occurrence of good days and 
poor days are equally likely to occur anywhere within the 
entire cycle then the distribution of these days should be 
uniform. Figure 3 illustrates this. A second statistical 
test should illustrate uniformity or non-uniformity. If the 
distribution is not uniform, this test should be able to 
identify the distributional form. Consequently, another 
statistical model should be used to illustrate uniformity 
or non-uniformity. 

C. QUALITY OF STATISTICAL MODEL 

When alternative statistical models are equally appro- 
priate, the technique with greater quality is chosen. Al- 
though there may be various ways to determine the quality 
of a model, the most objective measure of quality is power- 
efficiency. Suppose two tests, A and B are equally appro- 
priate. If test A requires samples to have the same power 
as test B has with samples then test B has 
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Frequency 




A = Theoretical envelope of biorhythmically predictable good 
days ; 

B = Hypothetical distribution of observed good days; 

C = Theoretical envelope of biorhythmically predictable poor 
days ; 

D = Hypothetical distribution of observed poor days. 

Figure 2. Theoretical Biorhythm Superimposed Upon Hypothetical 
Distribution of Observed Good and Poor Performance 
Days . 

/ N 

POWER-EFFICIENCY = (100) f 

where power is defined as the probability of rejecting the 
null hypothesis when it is false [Siegel, 1956]. 

Considering the characteristics of the data available, 
and appropriateness and quality of the statistical models, 
the Chi-Square Contingency Test for independence of criteria 
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A 



Frequency 



> 




time 



B 



A = Uniform distribution of observed good days; 

B = Uniform distribution of observed poor days. 

Figure 3. Hypothetical Sample from a Uniform Distribution 
of Good and Poor Performance Days. 

and the Chi-Square Goodness -of -Fit-Test for distributional 
form were chosen as primary statistical techniques. Secon- 
dary techniques are available to give supporting evidence. 
When sample size is so small that the Chi-Square Contingency 
Test's expected value criterion does not hold, the Fischer 
Test for exact probability (without Tocher's modification) 
should be used. When the Chi-Square Goodness -of -Fit -Test ' s 
expected value criterion does not hold and category combina- 
tion reduces the degrees of freedom to an undesirable level 
the Kolmorgorov- Smirnov Test for distributional form should 
be used. 
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These statistical techniques indicate the presence or 
lack of presence of a relationship and define the distribu- 
tional form. They do not attempt to explain why the rela- 
tionship is or is not present. Elaborate linear and 
curvilinear regression techniques are available for the de- 
velopment of quantitative explanatory models. 

In addition to analyzing the data with respect to de- 
pendency and distributional form, it would also be desirable 
to measure the degree of correlation between the theoretical 
biorhythm and the actual observations. This can be most 
easily accomplished by utilizing the Chi-Square Contingency 
Coefficient. 

D. CHI-SQUARE CONTINGENCY TEST 
1 . Theory 

A two-way Chi-Square Contingency Test was used to 
illustrate the dependency between the day of occurrence of 
the observed event and the sequential day of the cycle. 
Biorhythm theory states that an individual's capacity to 
perform is dependent upon the sequential day of the cycle. 
Specifically, there is a greater capacity for an individual 
to excel during the first half cycle and a reduction in 
capacity during the last half of the cycle. Therefore, the 
Chi-Square Contingency Test was used to investigate the de- 
pendency between the day of occurrence of observed human 
performance and the sequential day of the cycle. 

a. Data is arranged in the following tableau: 
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criteria 1 



1 

2 

3 

criteria 4 
2 5 



1 2 3 4 5 



x . . 



rc 



b. Given a sample of n observations. 

c. Each observation can be classified according to 
two different criteria. One criteria has c levels with 

c > 2 and the other criteria has r levels with r > 2. 



c 

Z 



r 

Z 



j = 1 i = 1 



x. . = n. 
ij 



e . c r 

Z Z P. . = 1. 
j=l i=l ^ 

f. nP. . > 5. 

il - 

g. Let x.. be the number of observations which fall 

il 

into class i criteria one and class j criteria two. 

h. 

falls into class i criteria one and class j criteria two. 



Let P^j be the probability that an observation 



l . 



P. 

l 



x . 

i ■ 

n 



where x. 



= Z x. . 
• n il 

1 = 1 J 



x . r 

— where x . = Z x. . 
n -I i=1 il 
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k. H : P. . = (P. ) (P .) for all i and i 
O 1J 1 • • J J 



H 1 : P. - f (P. ) (P -) 
1 ij v 1 • J K • j J 



1. 



If H is true: 
o 



c r 
V = E E 
j=l i=l 



C X ij - nP i. P .j)' 
nP- P . 

i- • j 



= X 2 (r-l)(c-l) 



(1-a) 



m. Reject H if: 
J o 



V > X 2 (r-l)(c-l) 



d-a) 



2 . Procedure for the Application to the Theory of 
Biorhythms 

a. Data tableau 



Criteria 1 



1 st is 
Cycle 



Criteria 
2 



2 nd % 
Cycle 



Observed Observed 

Good Days Bad Days 



X 11 


X 12 


X 2 1 


X 22 



s t 

1 *2 cycle: represents biorhy thmically predictable 

good performance days 

2 n< ^ h. cycle: represents biorhy thmical ly predictable 

poor performance days . 



are 

X 21 



b. x.. represents the number of observations which 
ij v 

categorized by both criteria levels i and j . For example 
represents the number of observed good performance days 
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I 






on which the individual's performance was above average, 
which .occurred during the second half of the cycle. 

Figure 4 is a graphical illustration of x^, x i2’ 



* 2 ^ , and x ^ . 

c. P. and P . represent the maximum likelihood esti- 

1 • *3 

mator for the probability that an observation is character- 
ized by the respective criteria level. For example, P^ t 
represents the probability that a good or bad day occurs 
during the 1st half of the cycle. 





2 






2 






X. 

1 • 


Z 

- 3 = 1 


x. . 
13 


P = 


Z x. . 
i-1 






n 


2 


2 


1. 

•J 


2 2 








Z 

3 = 1 


Z x. . 
i-1 1J 




Z Z x. . 

j=l i-1 ^ 








X 11 


+ V 

“12 


P - 


X 11 + 


X 2 1 




x n 


+ X 12 


+ x 2 1 + X 22 


^•1 


X 11 + X 12 + 


X 2 1 + 


X 22 




X 2 1 


+ X 2 2 


P = 


X 12 + 


x 22 




x n 


+ X 12 


+ X 21 + X 22 


•2 


X 11 + X 12 + 


x 2 1 + 


X 22 



e. For best results 

E(x. 0 = nP. P • > 5. 
^ i y i • • 3 - 

f. If H is true: 

o 



g- 



V = 



2 

Z 



2 (x. . 
Z — il 



nP. P .) 

i • • 3 



j=l i-1 



nP. P . 

i* *1 



Reject II if: 

V > X 2 (r-l)(c-l) (1 ' a) . 



X 2 (r-l)(c-l) (1_a) 
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Theoretical Biorhythm superimposed upon hypothetical dis- 
tributions of observed good performance and poor perform- 
ance days . 



1st h cycle 2nd h. cycle 




Xn = total number of observed good performance days occurring 
during the 1st % cycle. They appear in quadrant II. 

X 12 = total number of observed poor performance days occurring 
during the 1st ^ cycle. They appear in quadrant III. 

^21 = total number of observed good performance days occurring 
during the 2nd % cycle. They appear in quadrant I. 

^22 = num ber of observed poor performance days occurring 

during the 2nd cycle. They appear in quadrant IV. 



Figure 4. Location of Statistical Data Within a Biorhythm. 

E. CHI-SQUARE COEFFICIENT OF CONTINGENCY 

The contingency coefficient is especially useful in that 
it measures the degree or strength of the relationship which 
was detected using the Chi-Square Contingency Test for de- 
pendence of criteria. It is calculated as follows. 
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X 

X 2 + N 



C = Coefficient of Contingency. 

X 2 = Test statistic value (V) calculated from 
the Chi-Square Contingency Test (previously discussed). 

N = Total number of observations. 

F. CHI-SQUARE GOODNESS-OF-FIT TEST 

If the occurrence of a categorical event was equally 
likely to occur anywhere within the cycle then the distribu- 
tion of these events would be uniform. Consequently, valida- 
tion of the theory of biorhythms necessitates illustrating 
that the occurrence of categorical events is not uniformly 
distributed. Investigating distributional form is most 
easily accomplished by utilizing the Chi-Square Goodness-of- 
Fit Test. 

1. The data is organized in the following tableau. 



Interval 


Frequency 


P. 

1 


nP. 

l 


a < x < 


b 


Y i 


P 1 


ECVj) 


b < x < 


c 


Y 2 


P 2 


E( y 2 ) 


g < x < 


h 


Y k 


' P k 


E(Y k ) 



2. n: represents the sample size. 

3. a < x < b: represents an interval of time. 

4. Y^: represents the frequency of occurrence of the 

observed phenomena within the i z interval. 
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5. : represents the probability that an observed 

event .will occur, under the null hypothesis, during the 
associated interval. 

6. nP^ : represents the expected number of observations 

within the i^ 1 interval if the distribution of the frequency 
of occurrence of the phenomena conformed to the null 
hypothesis. 

7. For most significant results nP^ > 5. 

8. H : Hypothesized distributional form 

o 

: The distribution is other than as hypothesized. 

9. If H is true: 

o 



V 



k 

E 

i= 1 



< Y i 



k-r-1 = 
k = 



r = 



- nP.) 2 
1 J 

nP . 



x 2 (k-r-1) 



(1-a) 



degrees of freedom for the 

number of intervals 

number of estimators used 
sized distribution. 



test statistic 



in the hypothe- 



10 . Reject H q if : 

V > X 2 (k-r-1) ( 1_a ) . 

11. In an application to biorhythms, the initial test is 

for the presence of a uniform distribution. 

II : the distribution is uniform, 

o 

: the distribution is other than uniform 

p = 

*1 h - a 

P = d ~ c 

y 2 h - a 



etc . 
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Acceptance of the Contingency Test's null hypothesis 
indicates that the day of occurrence of observed events is 
independent of the sequential day of the cycle. In other 
words,- there is no relation between the days of occurrence 
of biorhythmically predictable events and the days of occur- 
rence of the observed events. This is an indication that 
the distribution of the occurrence of sampled events should 
be uniform. It is considered conclusive evidence (statis- 
tically) that the Theory of Biorhythms is invalid, if for 
the same data, H q for the Chi-Square Contingency Test and 
the Chi-Square Goodness -of -Fit Test is accepted. Rejection 
of both null hypotheses suggests validity of the theory. 
Further investigation of the distributional form would be 
desired. The presence of a known distributional form such 
as a normal or poisson would be a useful predictive device. 

G. ESTIMATION OF PARAMETERS 

Investigating the presence of more complex distributional 
forms usually results in estimating parameters of the hypothe- 
sized distribution. Consequently, estimators with good 
properties are important. The method of moments and the 
method of maximum likelihood techniques were considered. 

1 . Methods of Moments 

a. Given a population random variable X is dependent 
upon k unknown parameters 0^, ®2 , ‘ , ‘ , ®k‘ 

b. In general the first k population moments will 
depend upon these parameters (i.e., E (X ) = g (0^ , 0^ > • • • , OjJ ) • 
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c. Equate the first k sample moments , . . . ,M^ , 

to the first k population moments 



1 n k 

M, = - E (X.) , 
k n . t v 1 J 
i = l 



d. Solve for the desired parameter. 

2 . Method of Maximum Likelihood 

a. Develop a likelihood function 

n 

= ^ f jr ( X ) • 

n 

b. Since .E 1 f (x) is just another way of writing a 
sample's joint distribution function, the value of f(x) is 
actually a probability. 

c. The idea of the MLE technique is to find the 
value of the unknown parameters within the joint density 
which will maximize it. This value maximizes the probability 
of occurrence. The procedure is as follows: 



L C 0 ) = ^ f (x) 

i=l x 



K (0) = In (L (0) ) 

d (K (0) ) _ d(ln(L(0))j _ 0 
dO d0 U 

A 

solve for 0 and label it 0. 



In most cases 
of maximum likelihood 
those cases for which 
the properties of the 



the method of moments and the ..t^thod 
yield the same results. However, in 
the parameter estimators are different, 
estimators should determine which 
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estimator will be used. In many cases it is desirable to 
have unbiased, consistent, estimators w ith the least vari- 
ance and mean square error (MSE) . 

. 3 . Bias 

is an unbiased estimator for 0 iff EC©) = 0. The 
bias in is given by b^ = E(© ) - 0. 

4 . Consistency 

is a consistent estimator for 0 

iff limit EC©) = 0 
n ->■ 00 

limit V ((h) ) = 0. 

n co 

5. Cramer Rao Lower Bounds 



nE[Cg| log f x (x)) 2 ] 

mse = — C 1 + b . ' CpJ . ) 2 # 

nE[(gf log f x (x)) 2 ] 

Table II summarizes the combined statistical model 
suggested for testing the presence of a biorhythm. 

H. CRITICALITY 

Critical days are defined to be days on which an indi- 
vidual's potential capacity for performance is unstable. If 
an individual encounters a challenging situation on a bic- 
rhythmically critical day then it is possible that he will 
either have extreme negative potential or extreme positive 
potential. for example, it has been suggested that heart 
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Statistical Model 



Application 



1. x 2 Contingency Test 

(Use Fischer Test 
for small samples) 



Tests for the presence of a 
relationship between the 
different levels of perform 
ance and the days on which 
they occurred. 



2. x 2 Coefficient of 
Contingency 



Tests the strength of the 
relationship detected in 
step #1. 



3. x 2 Goodness -of -Fit Test 

(Use Kolmogornov- 
Smirnov Test for 
small samples) 



1st Application: 

Test for uniformity (test 
for the presence of cyclic 
phenomena) . 

2nd Application: 

Test for hypothesized dis- 
tributional form. This 
allows the analyst to in- 
vestigate the following: 

a) Form of cyclic phenomena 

b) Periodicity 

c) Frequency 

d) Criticalities 



Note : 

Tests should be applied in the above order. 

Detailed discussion and application procedures for the 
Fischer Exact probability Test and the Kolmogornov-Smirnov 
Test for distributional form can be found in Siegel [1956]. 



Table II. Combined Statistical Model for Testing the Presence 
of a Biorhythm. 
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attack victims encountered their attacks on critical days. 

It is obvious that a heart attack can be classified as a 
critical event. However, did a statistically significant 
number of heart attack victims experience their heart at- 
tacks on biorhythmic critical days? Thus, criticality be- 
comes an important factor when investigating the theory. 

For this reason, probabilistic arguments have been formu- 
lated for use in the following statistical model presented 
for investigating the significance of an event occurring on 
a critical day. 

1 . De finitions 

Single Criticality = day on which only one biorhythm 
crosses the time axis. 

Double Criticality = days on which any two biorhythms 
intersect while crossing the time axis. 

Triple Criticality = day on which all three biorhythms 
intersect while crossing the time axis. 

P(PCD) = probability that a physical critical day 
occurs by chance. 

P(ECD) = probability that an emotional critical day 
occurs by chance. 

P(ICD) = probability that an intellectual critical 
day occurs by chance. 



P(PCDxECD) = probability that a double critical day 
occurs by chance where the intersecting biorhythms are physi- 
cal and emotional. 
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P(PCDxICD) = probability that a double critical day 
occurs by chance where the intersecting biorhythms are physi- 
cal and intellectual. 

P(ECDxICD) = probability that a double critical day 
occurs by chance where the intersecting biorhythms are emo- 
tional and intellectual. 

P (PCDxECDxICD) = probability that a. triple critical 
day occurs by chance. 

It is possible that there are differing degrees of 
double criticalities. For example, it may be that a physi- 
cal x emotional double critical day possesses a higher 
degree of criticality than does a physical x intellectual 
double critical day. 

2 . Occurrence of Critical Days 

The following table illustrates the frequency of 
occurrence, in days of all types of critical days. 



Criticality Biorhythm 



Frequency (days) 



Single Physical (P) 11.50 
Single Emotional (E) 14.00 
Single Intellectual (I) 16.50 
Double P x E 161.00 
Double P x I 189.75 
Double E x I 231.00 
Triple PxExI 2656.50 



Table III. Frequency of Occurrence of Criticalities. 



The frequency (days) is based upon the respective 
cycle length. For example, a single criticality for a physi- 
cal biorhythm occurs every 11% days. Since the physical 
biorhythm is hypothesized to be a 23-day sinusoid, it crosses 
the axis every 11% days (i.e., 23/2 = 11%). 
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If the events are equally likely (uniform) to occur 
at any point in time, then the following probabilities mea- 
sure the likelihood that the observed event will occur on 
the given biorhythmic critical day. Calculation was based 
on the frequency of occurrence of critical days and con- 
sideration was given to the uncertainity of birth data. 



Chance Probability 


0 Hr 

Variance 


i 6 Hr 
Variance 


± 12 Hr 
Variance 


P (PCD) 


.0869 


. 1304 


. 1739 


P(ECD) 


.0714 


. 1071 


. 1428 


P(ICD) 


.0606 


.0909 


. 1212 


P (PCD x ECD) 


.0062 


.0093 


.0124 


P (PCD x ICD) 


.0053 


.0078 


.0105 


P (ECD x ICD) 


.0043 


.0065 


.0086 


P (PCD x ECD x ICD) 


.00037 


.0005 


. 0007 



Table IV. Chance Probability of Occurrence of Criticalities. 

3 . Selection of Statistical Model for Criticality 
Biorhythm data is naturally dichotomous into those 

days on which critical events occurred and those days on 
which critical events did not occur. Since an important 
question to be answered in validating the concept of criti- 
cality would be; "Are the number of observed critical days 
statistically significantly greater than the number that 
would be expected if critical event days were to have oc- 
curred by pure chance?"; the binomial model was chosen. 

4 . Binomial Model 

Assume X~Binomial (n,p) 

H : P = *P(critical day occurs by chance) 

H^: P > P (critical day occurs by chance) 

a : Level of significance 
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^Select appropriate probability from chance prob- 
ability table. 

Suppose that x observations out of n possible ob- 
servations are classified as critical events. 

Define P (X < x) = probability that at most x criti- 
cal events could have occurred by chance. 

n 

P(X < x) = E (?) P^l-P) 11 " 1 
i=0 1 

Decision: If P < a reject H in favor of H, . 

The decision to reject H q says that the number of 
critical events which have occurred is significant. 

For n > 25 it is suggested that a normal approxima- 
tion to the above probability be used. 

P(X < x) = 4> ( ^ x± ' 5) ~ nP ) . 

/nP(l-P) 

Additional information: 

E(X) = expected number of critical events under 

H = nP. 
o 

V(X) = variance of critical events under H q = nP(l-P). 

The Binomial Model is primarily used when investigat- 
ing the significance of single criticalities. The Multinomial 
Model should be used when investigating the statistical sig- 
nificance of data which contains single, double, and/or triple 
criticalities . 

5 . Multinomial Model 

Suppose that N observed events are partitioned into S 
mutually exclusive events Aj,A 2 ,...,A s with' respective 
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probabilities of occurrence P^,P 2 ,...,P S such, that + 

A- + . . . + A = N and P, + P- + . . . + P =1. Then P(A, 
2 s 12 s v 1 

occurs k.. times; A_ occurs k„ times; ...; A occurs k 
1 ’2 2 s s 

times) is equal to 





p^l p^2 p^3 

1 2 3 




I. BIORHYTHM ALGORITHM 



1. Note the date (month-day-year) of the observed event. 

2. Calculate the subject's age in days, from the sub- 
ject's birth date up to and including the date of the ob- 
served event. The calculation must include regular leap 
year and centurial leap year considerations. 

3. Let 

A = subject's age in days, 

CL = cycle length = 23 days for physical cycle 

28 days for emotional cycle 
33 days for intellectual cycle 

CC = whole number of completed cycles 

(whole number associated with A/CL) 

SD = the day of the current biorhythm cycle on which 
the observed event occurred. 



SD = ( cT - CC) ( CL ) * 

J. DATA COLLECTION 

Statistical validation of the theory necessitates per- 
forming two complete statistical analyses. The first analysis 
should use empirical data collected in a laboratory environ- 
ment, for each subcategory. Experimental design should screen 
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out as much exterior influences and as much interaction as 
possible, from the subcategorical data being collected. 

The second analysis should utilize actual field data. It 
is here that cycle interaction will have to be assessed. 

For example, if the analyst chooses to gather bowling data 
for physical field data, he should be aware that there 
exists some degree of emotional and intellectual interac- 
tion. If the analyst chooses to analyze academic test 
scores, he should be aware of which tests challenged the 
subject's intellect, and that the importance of the test and 
the length of the test (i.e., time) will cause emotional 
and physical interaction respectively. 

Since human behavior is extremely complex, there are a 
multitude of influencing factors. Therefore, it should be 
understood that biorhythms are not 'cure-alls'. A biorhythm 
is hypothesized to be merely one very small indicator of 
the variation of performance potential. The many influenc- 
ing factors can easily contaminate biorhythmic analysis. 
Consequently, it is not expected that field data will yield 
as strong a statistical significance as experimental data. 
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V. ANALYSIS OF DATA 



A. ANALYSIS OF INTELLECTUAL DATA 

Four graduate students at the Naval Postgraduate School 
were chosen as subjects for investigation of intellectual 
biorhythms. Academic grades from 15 courses measured over 
a 14 month period provided 115 data points for analysis. 

Each student, pursuing his Master's Degree in Operations 
Research, was evaluated with respect to the same courses. 

To eliminate a degree of subjectivity, only mathematically 
oriented courses were used. It was believed that grades 
were arrived at more objectively and that the subject's 
identification of his own neutral, average performance 
would be more accurate. Identification of this average per- 
formance is considered critical since it allows the analyst 
to categorize those performances which were above and those 
performances which were below average. Table V summarizes 
the data base. 

The data was next arranged in the Intellectual Data 
Tableau which is shown as Table VI. 

1 . Tests for the Presence of an Intellectual Biorhythm 
a. x 2 Contingency Test 

II : Observed good performance days and observed 

poor performance days are independent of where they occur 
within the intellectual biorhythm. 

: Observed good performance days and observed 

poor performance days are dependent of where they occur with 

the intellectual biorhythm. 
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Ranking of 

Day of Current Performane 

Code Date of Exam Grade Ach/Grade Poss Intellect Cycle Potential 
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Table V. Intellectual Biorhythm Data Base. 



Subject #1 Table V. continued Ranking of 

Day of Current Performance 
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Subject If 3 Ranking of 

Day of Current Performance 

Course Code Date of Exam Grade Ach/Grade Poss Intellect Cycle Potential 

MC 12 Oct 72 72/100 00 
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Subject #4 Ranking of 

Day of Current Performano 

Course Code Date of Exam Grade Ach/Grade Poss Intellect Cycle Potential 
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Table V. continued 

Sub j cct ,4 Ranking of 

Day o£ Current Performance 

Course Code Date of Exam Grade Ach/Grade Poss Intellect Cycle Potential 
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cuius exam. If the class average was 75 and the subjects overall average 
exam was 80, then he was ranked ++. Since many ranking scales used in the 
are ordinal, at best, the choice between + and ++ is subjective. 



Avg Slightly 
Cycle Per Above 
Day 0 Avg (+) 



Well 
Above 
Avg (++) 



Slightly Well 
Below Below 

Avg (-) Avg (--) 



Critical 0 



1 1 

2 13 

3 1 

4 

5 3 

6 2 

7 11 

8 3 3 

9 2 

10 1 

11 4 

12 

13 2 

14 1 

15 1 1 

16 

17 

18 1 1 

19 

20 13 

21 

22 1 

23 2 1 

24 1 

25 2 

26 2 

27 1 

28 11 

29 2 

30 1 1 

31 2 

32 1 

33 



1 

1 

2 

1 

1 

1 



4 

1 

1 

1 



1 



2 



1 



1 

2 

2 

1 



2 



1 

2 

3 

3 

2 

2 

1 

2 

3 

1 

1 

2 

1 

2 



1 



1 



1 



1 

1 

1 



1 



1. The data points within the tableau represent the frequency 
of occurrence. For example, there were four above average 
performances on the 11th day of the cycle. 

2. The data can now be easily categorized for use in the 
various statistical models. 



Table VI. Intellectual Data Tableau. 
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a : Level of Significance = .005. 





Observed Good 
Performance 
Days 


Observed Poor 
Performance 
Days 


1st % Cycle 


21 


4 


2nd k Cycle 


15 


9 


3rd % Cycle 


8 


14 


4th k Cycle 


2 


16 



Xj ( . 995) = 12.8 V = 25.7 



Table VII. Intellectual Data x 2 Contingency Test. 

Decision: Reject H in favor of IL . This in- 

dicates the presence of a relationship between the different 
levels of performance and the biorhythmic days on which they 
occurred . 

b. x 2 Coefficient of Contingency 

^ + N \25.7 + 89 ,46 ' 

This represents the strength of the relation- 
ship. It is important to realize that C = .48 indicates 
that the relationship is strong. Notice that even if the 
relationship were perfect, it would be impossible for C = 1, 
since V/V + N will always be less than 1. 

2 . Distribution of Average Performance Days 
a. Test For Uniformity 

X 2 Goodness-of -Fit Test 

H : The distribution of average performance day 

is uniformly distributed over the intellectual biorhythm. 
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H^: The distribution of average performance 

days is not uniformly distributed over the intellectual 
biorhythm. 

a : Level of Significance = .01. 



Interval (days) Frequency nP^ 



1 


< 


t 


< 


8 


5 


8/33 


6.06 


8 


< 


t 


< 


16 


4 


8/33 


6.06 


16 


< 


t 


< 


26 


8 


10/33 


7.57 


26 


< 


t 


< 


33 


8 


7/33 


5.30 



x| ( . 99) = 11.3 V = 2.28. 



Table VIII. Intellectual Data (Average Performance Days): 
X 2 Goodness-of-Fit Test for Uniformity. 

Decision: Do not reject H . 

J o 

Kolmogorov-Smirnov Test 

Utilizing the same H q , , and a as above: 



Interval 


1 < t < 8 


8 < t < 16 


16 < t < 26 


26 < t < 33 


Frequency 


5 


4 


8 


8 


V x) 


8/33 


16/33 


26/33 


33/33 


S 25 (x) 


5/25 


9/25 


17/25 


25/25 


l F o ' S 2 5 1 


.0424 


. 1248 


. 1079 


00.00 



Critical D = .32 Calculated D = .1248. 



Table IX. Intellectual Data (Average Performance Days) : 
Kolmogorov-Smirnov Test for Uniformity. 

Decision: Do not reject H . 

J o 

Since the distribution is uniformly distributed, 
there is no need to further investigate distributional form. 
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The above analysis confirms the assumption that average per- 
formance days can occur anywhere within the cycle with equal 
probability (uniformity) . 

3 . Distribution of Above Average Performance Days 

For stronger test results, the slightly above aver- 
age (+) and well above average (++) performance days were 
combined for analysis under the heading of above average 
performance days. 

a. Test for Uniformity 

X 2 Goodness -of- Fit Test 

H q : The distribution of above average perform- 

ance days is uniformly distributed over the intellectual 
biorhythm. 

H^: The distribution of above average perform- 

ance days is not uniformly distributed over the intellectual 
biorhythm . 



a : Level of Significance = .01. 

Interval (days) Frequency nP^ 



1 


< 


t 


< 5 


13 


5/33 


6.96 


5 


< 


t 


< 10 


10 


5/33 


6.96 


10 


< 


t 


< 15 


13 


5/33 


6.96 


15 


< 


t 


< 20 


4 


5/33 


6.96 


20 


< 


t 


< 25 


4 


5/33 


6.96 


25 


< 


t 


< 33 


2 


8/33 


11.20 








Xc ( • 99) 


= 15.1. 


V = 21.88. 





Table X. Intellectual Data (Above Average Performance Days) 
X 2 Goodness- of -Fit Test for Uniformity. 



Decision: Reject H in favor of H 
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Kolmogorov-Smirnov Test 



Utilizing 

reviously : 


the same 


H , H 1 , 
o’ 1 J 


and a as 


discussed 


Interval 


l<t<5 


5<t<10 


10< t<l 5 


15<t< 20 


20<t<2 5 


2 5<t< 33 


Frequency 


13 


10 


13 


4 


4 


2 


F 0 « 


5/33 


10/33 


15/33 


20/33 


25/33 


33/33 


s 46 (>0 


13/46 


23/46 


36/46 


40/46 


44/46 


46/46 


|F - S| 


.1311 


.1970 


.3281 


.2635 


. 3504 


00.00 




Critical 


D = . 2403 Calculated D 


= .3504 





•Table XI. Intellectual Data (Above Average Performance Days): 
Kolmogorov-Smirnov Test for Uniformity. 

Decision: Reject H q in favor of H^. 

Both of the above tests tested for uniformity. 

Lack of uniformity suggests the presence of a cyclic phenomena. 

b. Test for Distributional Form 

X 2 Goodness-of -Fit Test 

H : The distributional form of above average per- 

0 

formance days is N(8.5,5). 

H.. : H is false . 

1 o 

a : Level of Significance = .01. 



Interval 


(days) 


Frequency 


P. 

l 


nP. 

l 


- CO < t 


< 3 


7 


. 1357 


6.24 


3 < t 


< 12 


24 


.6226 


28.64 


12 < t 


< CO 


15 


. 2417 


11.12 



Table XII. Intellectual Data (Overall Above Average Perform- 
ance Days): x 2 Goodness -of -Fit Test for Distri- 

butional Form. 
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V = 2.19 



X^(.99) = 9.21 

Decision: Do not reject H . 

J o 

4 . Distribution of Below Average Performance Days 

For stronger test results the slightly below aver- 
age (-) and well below average (--) performance days were 
combined for analysis under the heading of below average 
performance days. 

a. Test for Uniformity 

X 2 Goodness-of -Fit Test 

H : The distribution of below average perform- 

ance days is uniformly distributed over the intellectual 
biorhythm. 

H^: The distribution of below average perform- 

ance days is not uniformly distributed over the intellectual 
biorhythm. 



a : Level of Significance = .10. 

• • 

Interval (days) Frequency nP^ 




Table XIII. Intellectual Data (Below Average Performance Days) 
X 2 Goodness-of-Fit. Test for Uniformity. 

Decision: Reject II in favor of FL . 

J o 1 
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Kolmogorov-Smirnov Test 

Utilizing the same H q , H^, and a as previously 

discussed 



Interval 


1 < t < 8 


8 < t < 


16 16 < t < 27 


27 < t < 33 


Frequency 


4 


9 


21 


9 


V- x > 


8/33 


16/33 


27/33 


33/33 


S 43 (x) 


4/43 


13/43 


34/43 


43/43 


V S 43 


. 1494 


. 1825 


.0275 


00.00 




Critical D 


= .1753 


Calculated D = 


. 1825. 



Table XIV,. Intellectual Data (Below Average Performance Days) : 
Kolmogorov-Smirnov Test for Uniformity. 

Decision: Reject H q . 

b. Test for Distributional Form 

X 2 Goodness -of-Fit Test 

H : The distributional form of below average 

performance days is 3(4,2) 

H 1 :H is false. 

1 o 

a : Level of Significance = .01. 

Interval (days) Frequency nP^ 



1 


< 


t 


< 


20 


13 


. 350 


15.050 


20 


< 


t 


< 


30 


24 


.575 


24.725 


30 


< 


t 


< 


33 


6 


.075 


3.225 



X2C.99) = 9.21 .V = 2.69. 



Table XV. Intellectual Data (Overall Below Average Perform- 
ance Days): x Z Goodness -of -Fit Test for Distri- 

butional Form. 

Decision: Do not reject H . 

J o 



55 



B. ANALYSIS OF ACCIDENT DATA 



Sixty-six accident claims from the Canadian Forest Prod- 
ucts Ltd. (Howe Sound Pulp Division) were analyzed primarily 
with respect to their physical biorhythms. Those accidents 
analyzed were single victim accidents in which the accident 
was not believed to be caused or influenced by a second 
party. Table XVI summarizes the data base. 

The data was next arranged in the Physical Biorhythm 
Data Base tableau which is shown as Table XVII. 

1 . Tests for the Presence of a Physical Biorhythm 
a. First Analysis 

Since accidents are considered in the below 
average performance category, the y 2 Contingency Test and 
X 2 Coefficient of Contingency are inappropriate. Conse- 
quently, the analysis begins with a test for uniformity. 

X 2 Goodness-of -Fit Test 

H : The occxirrence of accidents is uniformly 

o } 

distributed over the physical biorhythm. 

H^: The occurrence of accidents is not uniformly 

distributed over the physical biorhythm. 

a : Level of Significance = .1. 

Interval (days) Frequency P. nP. 



1 


< 


t 


< 


5 


10 


5/23 


14 . 35 


5 


< 


t 


< 


10 


08 


5/23 


14.35 


10 


< 


t 


< 


15 


20 


5/23 


14.35 


15 


< 


t 


< 


20 


15 


5/23 


14.35 


20 


< 


t 


< 


23 


13 


3/23 


8.60 



Xj(.90) = 7.78 V = 8.63 

lable XVIII. Accident Data (1st Analysis): x 2 Goodness-of- 
Fit Test for Uniformity. 
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bj< 

# 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



Day of Current 
Phys Biorhythm 



21 

19 

13 

11 

03 
10 
21 
00 
19 
08 

19 
22 

20 

13 
07 
02 

04 

14 
16 
22 
16 

10 
1 1 

05 
05 
00 
14 
16 
21 
18 

19 
10 

20 
03 
22 
12 
00 

13 
21 
12 
17 

14 
21 
07 
11 



Day of Current 
Emot Biorhythm 



22 

19 
22 
02 
02 
22 
21 

20 
11 
25 
11 
13 
05 
02 
11 
02 
24 

23 

24 
01 
15 

17 

18 
13 

09 
20 
02 
27 
27 

10 

09 
20 

19 
04 

10 

25 
01 
04 

20 
00 
02 
02 
10 
09 
19 



Day of Current 
Int Biorhythm 



24 
26 

25 
18 

09 
14 
19 

17 
19 

19 
04 

18 

23 
21 

20 

10 

24 
21 
16 
07 

04 

07 
09 
21 

05 
17 
05 

25 
19 
24 

08 

05 

30 
14 

07 

06 
01 
21 

31 
27 
31 
11 

26 

08 
26 



Table XVI. Accident Data Base. 
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Table XVI. continued 



Subject 

# 


Day of Current 
Phys Biorhythm 


Day of Current 
Bmot Biorhythm 


Day 

lnt 


of 

Bi 


46 


20 


11 




01 


47 


14 


17 




24 


48 


02 


03 




24 


49 


21 


09 




17 


50 


19 


11 




32 


51 


01 


15 




16 


52 


15 


10 




28 


53 


17 


09 




04 


54 


15 


23 




04 


55 


00 


18 




24 


56 


11 


15 




16 


57 


13 


01 




07 


58 


15 


11 




17 


59 


17 


26 




06 


60 


09 


14 




15 


61 


07 


15 




12 


62 


14 


25 




15 


63 


01 


16 




07 


64 


01 


09 




04 


65 


13 


18 




51 


66 


14 ' 


04 




20 



Current 

orhythm 
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Cycle Day 



Accident Frequency 



1 3 

2 2 

3 2 

4 1 

5 2 

6 

7 3 

8 1 

9 1 

10 3 

11 4 

12 2 

13 5 

14 6 

15 3 

16 3 

17 3 

18 1 

19 5 

20 3 

21 6 

22 3 

00 4 



Frequency refers to the total number of accidents which 
occurred on the respective biorhythm day. For example, of 
the 66 subjects, three experienced accidents on the first 
day of their physical biorhythm, two experienced their ac- 
cidents on the second day, etc. The data can now be easily 
categorized for use in the various statistical models. 



Table XVII. Physical Biorhythm Data Base. 
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Decision: Reject H in favor of FL . 

o 1 

Kolmogorov-Smirnov Test 

Utilizing the same , H^, and a as previously 

discussed 



Interval (days) 


i<t<5 


5<t<10 


10<t<15 


15<t<20 


20<t<23 


Frequency 


10 


8 


20 


15 


13 


F 0 (x) 


5/23 


10/23 


15/23 


-20/23 


23/23 


S 66« 


10/66 


18/66 


38/66 


53/66 


66/66 


V S 66 


.0655 


.1620 


.0764 


.0665 


00.00 



Critical D = .1502 Calculated D = .1620. 

Table XIX. Accident Data (1st Analysis) : Kolmogorov-Smirnov 

Test for Uniformity. 

Decision: Reject H in favor of H. . 

J o 1 

Rejection of the uniformity hypothesis by both 
statistical tests suggests the presence of a cyclic phenomena, 
b. Second Analysis 

The Theory of Biorhythms has a built in variance 
of ±12 hours (24 hours) due to the uncertainty of the exact 
hour of birth. Consequently, it is possible that the original 
calculations may be as much as 24 hours away from where they 
appear in the original data tableau. This analysis investi- 
gates the statistical significance in shifting the data 24 
hours (optimally) . 

X 2 Goodness -of -Fit Test 

H : The occurrence of accidents is uniformly 

o 

distributed over the physical biorhythm. 
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Cycle Day 



Accident Frequency 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 



1 

1 

2 

1 

1 

1 

2 

2 

1 

1 

2 

6 

5 

6 
3 
3 
3 
1 

5 
3 

6 
3 
7 



Table XX. Optimal Data Tableau. 



H^: The occurrence of accidents is not uniformly 

distributed over the physical biorhythm. 

a : Level of Significance = .005. 



Interval (days) Frequency nP^ 



1 


< t 


< 


5 


6 


5/23 


14.35 


5 


< t 


< 


10 


7 


5/23 


14.35 


10 


< t 


< 


15 


22 


5/23 


14.35 


15 


< t 


< 


20 


15 


5/23 


14.35 


20 


< t 


< 


23 


16 


3/23 


8.60 




xic. 


995) = 14.9 




V = 19.71. 





Table XXI. Accident Data (2nd Analysis): x 2 Goodness-of- 

Fit Test for Uniformity. 
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Decision: Reject H in favor of H., . 

J o 1 

Kolmogorov-Smirnov Test 

Utilizing the same H , and a as previously 

discussed 



Interval 


l<t<5 


5<t<10 


10<t<15 


15<t<20 


20<t<23 


Frequency 


6 


7 


22 


15 


16 


F „w 


5/23 


10/23 


15/23 


20/23 


23/23 


S 66« 


6/66 


13/66 


35/66 


50/66 


66/66 


F o ' S 66 


.1265 


.2378 


.1219 


. 1120 


00 . 00 


Critical D 


= .2006 


Calculated D = 


.2378 



Table XXII. Accident Data (2nd Analysis);' Kolmogorov-Smirnov 
Test for Uniformity. 

Decision: Reject H . 

The above analysis shows that shifting the oc- 
currence within the cycle by 24 hours yields more significant 
results. However, knowledge of the subject's exact hour of 
birth eliminates the need for shifting the data 24 hours and 
yields more reliable results, 
c. Third Analysis 

This analysis dichotomizes the data into loss- 
time accidents and non- loss - time accidents. A loss- time 
accident is an accident which results in the employee re- 
maining out of work for at least one day. The following 
tableau contains data from loss - time - accident subjects. 
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Sub j ect 


Day of Current 


Day of Current 


Day 


of Current 


# 


Phys Biorhythm 


Emot Biorhythm 


Int 


Biorhythm 


01 


*21 


22 




24 


11 


19 


11 




04 


18 


14 


23 




21 


20 


*22 


*01 




07 


21 


16 


*15 




04 


23 


*11 


18 




09 


39 


*21 


20 




*31 


41 


17 


*02 




*31 


53 


17 


09 




04 


58 


15 


11 


■ 


*17 


66 


14 


04 




20 



*It is possible, assuming ±12 hours variance, 
that these days could have been critical days. 



Table XXIII. Loss-Time Accidents. 



X 2 Contingency 

H : Loss-time 

o 

dents independent of where 
biorhythm. 



Test 

accidents and non- loss - time acci- 
they occur within the physical 



H^: Loss-time accidents and non- los s -time ac- 

cidents are dependent of where they occur within the physical 
biorhythm. 



a : Level of Significance = .05. 







Loss -Time 


Non-Loss -Time 






Accidents 


Accidents 


1st 


% Cycle 


1 


21 


2nd 


% Cycle 


11 


33 



X 2 (.95) = 3.84 V = 4.18. 



Table XXIV. x 2 Contingency Test. 

Decision: Reject H in favor of H 

J o 



This result illustrates more severe accidents 



are dependent upon the sequential day of the subjects 
biorhythm . 



X 2 Coefficient of Contingency 




V 



V + N 




C 



.244. 



This represents the degree of association or the 



strength of the relationship previously measured by the y 2 
Contingency Test. The Coefficient of Contingency is just 
another way of examining the strength of the rejection in 
the x 2 Contingency Test. Had the x 2 Contingency's null 
hypothesis (that of independence) has been rejected, there 
would have been no reason to calculate a Contingency Coeffi- 
cient since the y 2 deviation statistic would have been so 
small that the resulting Coefficient would have been insig- 
nificant. This can be illustrated by examining the formula. 

In order not to reject the null hypothesis, the deviation 
statistic would have to be less than 3.84. Assume that this 
statistic was 1.5. With the same sample size of 66, the 
Coefficient of Contingency is equal to .149 which is insig- 
nificant at any alpha level. A close examination of the 
parameters V and N illustrates that one must be careful when 
interpreting this coefficient. A plot of the Contingency 
Coefficient against sample size, holding the deviation constant 
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will show that as the sample size becomes large, the coeffi- 
cient becomes asymptotic to zero. Consequently, it is more 
sensitive to smaller sample sizes, as is the case with this 
biorhythm data. 

2 . Criticality 

Lack of data eliminates an indepth analysis of 
criticality, however, it is interesting to observe the 
occurrence of criticalities among loss- time accident victims. 

Let 

X = criticalities 
P (PCD) = .1739 
P(ECD) = .1428 
P(ICD) = .1212. 



Biorhythm 


E (X) 


V(X) 


# Critical Observations 


Physical 


1.91 


1.58 


4 


Emotional 


1.57 


1.35 


3 


Intellectual 


1.33 


1. 16 


3 



Table XXV. Accident Data Criticality Tableau. 

The following probabilities were calculated using 
the binomial model. 

Let 

X = # criticalities observed within a physical 
biorhythm out of 11 possible observations. 

P (PCD) = .1739. 

Table XXVI illustrates that the expected (mean) 
number of observed criticalities (1.91) occurs with highest 
probability while four or more criticalities can be considered 
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a rare event. Figure 5 illustrates that 90% of the proba- 
bility of occurrence of critical observations is contained 
within 1.5 standard deviations of the mean. 



n P (X=n) 

0 .123300000 

1 .283109200 

2 .297651200 

3 .187889600 

4 .077992200 

5 .023344900 

6 .004905300 

7 .000722300 

8 .000074400 

9 .000003700 

10 .000000158 

11 .000000003 



Table XXVI. Criticality Prob Tableau. 
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Figure 5. Observed Accident Data Density Function. 
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VI. DISCUSSION 



Results based on statistical models are heavily dependent 
on the researcher’s interpretation of the theory. One who 
considers criticality as the backbone of the theory will 
stress different aspects of biorhythmic analysis than will 
a researcher who believes criticality plays a minor role in 
biorhythmic prediction. Consequently, the choice and de- 
velopment of the combined statistical model illustrates that 
every aspect of a biorhythm should be analyzed prior to draw- 
ing conclusions. Use of the models in a thorough and rigid 
analysis would allow one to draw objective conclusions based 
on fact instead of formulating broad generalizations based 
on subjective interpretations of the theory. 

The analysis of intellectual and accident data are not 
meant to validate the theory. The analyses both illustrate 
the use of the statistical models as well as support this 
paper's interpretation of biorhythms. The analysis of in- 
tellectual data rigidly adhered to the combined statistical 
model for testing the presence of a biorhythm. At the .005 
level of significance the Chi-Square Contingency Test es- 
tablished that a relationship existed between the different 
levels of performance and the days on which they occurred. 
Supporting this result was an impressive .48 coefficient of 
contingency. As was expected the application of both the 
Chi-Square Goodness -of -Fit Test and the Kolmogorov-Smirnov 
Test (KS) to average performance data showed that days on 
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which the subject was not intellectually challenged (average 
performance days) were uniformly distributed over the entire 
cycle. This result is coincidental with biorhythm theory. 
Biorhythms measure variation in performance potential. On 
average days variation is zero and thus unpredictable and 
completely random as in the case of the uniform distribution. 
However, repeated application of the Chi-Square and KS Tests 
to the above and below average performance days illustrated 
non -uni f ormity which, according to the model implies a 
cyclic phenomena. Identification of the distribution of 
both the above and below average performance days was pos- 
sible. Identifying the distribution is instrumental in that 
it allows one to determine expected frequencies and their 
variations. Although application of both statistical tests 
was not necessary, the use of the additional test serves as 
a check and lends supporting evidence to the original analy- 
sis. For small samples the KS Test is more powerful since 
the Chi-Square Test is forced to combine categories when 
expected value criteria are not met. It is realized that 
the limited sample size presents difficulties when attempt- 
ing to dichotomize in attempts to gain as much information 
as possible from the data. However the significance of the 
results is impressive regardless of the sample size. The 
presence of only one critical event prevented a criticality 
analys is . 

There are many problems encountered when attempting to 
analyze data such as academic grades. It has been the 
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author's experience that exams on the graduate level are 
designed to challenge the student. Most conscienious stu- 
dents over-prepare for exams thus reducing the exams’ chal- 
lenge. This makes it difficult to collect and analyze clean 
data. If the intellectual biorhythm actually measures 
memory potential then a knowledge of potential on exam days 
may be a determinant in how much preparation is. required 
for the exam. 

In general, analyses of accident data is not considered 
appropriate in analyzing biorhythms since accident data 
isolates the analysis to the low performance potential area 
of the biorhythm. This single class of data prevents the 
use of the entire statistical model. For example the Chi- 
Squarc Contingency Test and the Coefficient of Contingency 
are not applicable as suggested. However application of 
both the Chi-Square Goodness-of-Fit Test and the KS Test 
established the existence of a non-uniform distribution of 
below average performance days. One may wonder what is the 
effect of considering a natural built in variance due to 
the uncertainty of the exact hour of birth and possibly the 
exact time of the observed event. Between critical days 
(above and below average performance) the variance has a 
negligible effect on whether the biorhythm is present. How- 
ever it has a profound affect on the criticality in that the 
assumption of variance increases the chance probability of 
an event occurring on a critical day. Further dichotomizing 
the time-loss and non- time- loss accidents allowed the assessment 
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of the severity of the accidents. The refinement in classi- 
fication resulted in applying the Chi-Square Contingency 
Test and the Coefficient of Contingency. This showed that 
time-loss accidents (more severe) were dependent upon the 
biorhythmic day on which they occurred. It is also sug- 
gested that the occurrence of critical events is signifi- 
cant. For a sample size of 11 the expected number of events 
occurring on a critical day is 1.91 with a variance of 1.58 
and a standard deviation of 1.25. The fact that four criti- 
cal observations, lie outside 1.5 standard deviations of 
the mean suggests the presence of the criticality phenomena. 
Although it was possible to calculate the exact probability 
of this event, it is felt that further attempts at accuracy 
will not be meaningful. 

P(3 phy critical days; 1 double critical day; 

7 non-critical phy days) = 

11 3 7 

_ n 7 C- 1739) ^C-0124) C. 8137) 7 = .02. 

As illustrated, this can be considered a rare event. On the 
average it is expected that 1.91 events would occur on criti 
cal days. It is also expected that, under a null hypothesis 
of pure chance, the expected number is not rigid and has a 
variance of 1.58 and a standard deviation of 1.25. This 
means that if the occurrence of critical events was pure 
chance one could expect as few as zero observed critical 
events or as many as 3.82. The fact that four critical 
events were observed suggests that criticality does not 
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occur in a pure chance fashion since more critical events 
occurred than were expected. Consequently it is felt that 
a criticality phenomena is present. 

Although both analyses do not comprise a validation, 
each analysis completely supports the theory in a very 
strong sense. 
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VII. CONCLUSION 



Although biorhythmic phenomena is considered esoteric 
and attributed to the occult sciences, the possibility that 
such isochronic phenomena exists should not be dismissed. 

As suggested there are scientific methodological ways of 
investigating the theory. The method presented here is not 
considered the only method. Preliminary investigations 
suggest the validity of the theory, however, only a scien- 
tifically based validation will be accepted by a scientific 
community . 
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